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Non-relativisitic second-order corrections to tlie wave function at origin in muonic and exotic 
atoms are considered. Tlie corrections are due to tlie electronic vacuum polarization. Such correc- 
tions are of interest due to various effective approaches, which take into account QED and hadronic 
effects. The wave function at origin plays a key role in the calculation of the pionium lifetime, 
various finite nuclear size effects and the hyperfine splitting. The results are obtained for the is 
and 2s states in pionic and muonic hydrogen and deuterium and in pionium, a bound system of ■k'^ 
and 7r~. Applications to the hyperfine structure and the Lamb shift in muonic hydrogen are also 
considered. 



INTRODUCTION 



electron mass Wg and in particular^ 



A number of atomic effects and, in particular, in ex- 
otic atoms involve in the non-relativistic approximation 
various local operators, proportional to the (5-function in 
coordinate space. The related contributions are propor- 
tional to the squared value of the wave function at origin 
|^nr(0)P. Two examples of such operators are the op- 
erator of interaction of the muon spin and the nuclear 
spin in a muonic atom (that is responsible for the hyper- 
fine structure) and the 7r"'"7r~7r°7r'^ vertex operator (that 
is responsible for the lifetime of the pionium atom). 

That is a common feature of various non-relativistic 
approximations and various effective non-relativistic ap- 
proaches, which are based on a separation of low-energy 
and high-energy physics. Physics of atomic scale con- 
tributes to the non-relativistic wave functions, while the 
higher energies and momenta are responsible for various 
contact terms. That is very much similar to the operator 
approach in theory of strong interactions. 

Meanwhile, there is an important difference between 
'conventional' atoms and various exotic atoms in the cal- 
culation of |^nr(0)P. In conventional (electronic) atoms 
the non-relativistic wave function is in most problems de- 
termined by its pure Coulomb value |\E'c(0)P and most of 
the corrections are either relativistic or have many-body 
origin. In contrast, in muonic and pionic atoms, there is 
a specific class of non-relativistic corrections, which can 
be still described by a non-relativistic potential. The or- 
biting mass TO in such atoms is much higher than the 



TO > TO^ ~ 207 TOe . 

The characteristic momentum in such atoms Zamc is 
higher or comparable to mgC, and thus the electronic 
vacuum polarization produces a non-relativistic potential 
with the radius of ^ h/{mec), which is somewhat larger 
than the atomic Bohr radius ^ h/{Zamc). 
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FIG. 1: Vacuum polarization corrections to electrostatic 
Coulomb potential (Vvp): the Uehling potential (V^p), the 



reducible two-loop (Vyp^'' 



(2)^ 



and irreducible {Vyp 
potentials (the Kallen-Sabry potential). 



two-loop 
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^ Technically, in the non-relativistic case that is the reduced 
mass that enters the equations. It may be somewhat below 
the muon mass rrtj,. The smallest values are in systems of 
Jlfi, TTfi and tt+tt": m^^ = 0.5 m^, m,r/j = 0.569... and 
nijrir = 0.660 . . . m„. 
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The related potentials depicted in Fig. [T] modify the 
value of the non-relativistic wave function at origin 

|*c(0)P^|*nr(0)|' = |*c(0)P 
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The related diagrams are presented in Fig. [2l The 
coefficient c\ for the one-loop corrections is due to the 
Uehling potential and was evaluated for a number of 
problems (see, e.g., [H, 0, 1 i, i, Q ) . 




FIG. 2: The first-order correction to the wave function at 
origin. The filled circle is for the (5(r). The double line stands 
for the non-relativistic reduced Coulomb Green function. 

The second-order effects (see Fig. ^ are due to subse- 
quent iterations of the Uehling term and due to Kallen- 
Sabry potential and here we present a calculation of C2 
for the Is and 2s states in muonic and pionic hydrogen 
and deuterium and for the pionium atom. 
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includes the Coulomb potential and its modification by 
the first- and the second-order vacuum polarization (see 
Fig. [J) 

Wp(r) = T/«(r) + V^^^\r) + V^^^{r) . 

We are to find the wave function ^'v(r) related to a po- 
tential l^(r'), or rather, only its value at origin. For the 
latter, we can introduce an additional perturbation 

and the obvious result of the perturbation theory, linear 
in (5V5(r), is 

^i;^ = a|*nr(o)P . 

That means that instead of finding the wave function we 
can calculate the perturbation theory expansion for en- 
ergy with the perturbation Vvp(?') + Vsij), taking terms 
linear in A and of a proper order in a. 

So, to find the coefficient c\ it is enough to consider 
linear terms in Vvp(?'), and for the latter to apply the 
Uehling term only. To find ci we have to evaluate two 
kinds of contributions, which are 

- terms quadratic in Vvp(»'), but including only the 

Uehling potential (Vvp(?') ^vp(?-)) (see Fig. Eb 
and d); 

- terms linear in Vvp(r), which include the second- 
order vacuum polarization (both reducible and irre- 
ducible: (Wp(r) -> yYV^\r) + vi§{r)) (seeFig.Hi 
and 6). 

The results of the calculation of the coefficients c\ and 
C2 are collected in Table HI 
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FIG. 3: The second-order vacuum polarization correction to 
the wave function at origin. 



II. GENERAL EXPRESSION FOR THE 
CORRECTION TO THE WAVE FUNCTION AT 
ORIGIN: ORDERS a AND 

We consider a non-relativistic two-body system, which 
interaction 

T/(r) = Vb(r)+Wp(r) (2) 



TABLE I: The results of the calculation of the coefficients 
Ci and C2 in various atoms for the Is and 2s states. The Ci 
coefficient was discussed in literature (see, e.g., P, 0, l^i 0, 
0! while the results for C2 are found in this work. The 
accuracy of the presented results is determined by rounding 
of the numerical results. 



Atom 


mr/rrie 


Is 


2s 


Cl 


C2 


Cl 


C2 




136.566 


1.35025 


5.4378 


1.13440 


4.3723 




185.841 


1.73115 


7.2558 


1.40425 


5.5552 




195.742 


1.80116 


7.6038 


1.45230 


5.7730 




237.764 


2.07748 


9.0209 


1.63850 


6.6402 


nD 


254.215 


2.17742 


9.5504 


1.70477 


6.9584 



More details of the calculation of the C2 coefficient can 
be found in Table [TTl where we present separately all 
contributions for the Is state. The result for the Fig. 
contribution is split into two terms. That reflects the 
fact that in general a contribution in the third order of a 
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perturbative theory is determined by the expression (see, 
e.g., 0,i) 



5V - AMI) 



(3) 



where AEnJ — {'i'ns\SV\"^ns}, SV is a sum of all pertur- 
bations under consideration and '^ns and G are the wave 
function and the reduced Green function, respectively, 
of the unperturbed problem (i.e., of the non-relativistic 
Coulomb problem in our case). 

A calculation of the subtraction term in ([3]), which is 
of the form 



A^^^ns) = -A^W X {^n,\5VG^SV\-i'r. 



(4) 



is different from that of the main term in ^ and in fact 
for the wave function it is even somewhat more compli- 
cated. 



TABLE II: Contributions to the value of C2(ls) for different 
diagrams in Fig. |31 The result for the Fig. [31; and d contribu- 
tions is split into two parts (of. (O). 



(cf. Fig. |4]) for the VP contributions, which is (see [lE 
for detail) 



AE^^ils) = (^2.61419^ + 12.54584 (^) j AeI^M , 

(6) 



A^^^(2s) = (2.3U5l^ + 10.65790 (^)'^ A£;|^"](2, 



(7) 



where 



AEl^i{ns)^-il + a,) 



a{Za) mc m ji /m^NS 



y 



(8) 

/iAT is the nuclear magneton, stands for the proton mag- 
netic moment, and is the reduced mass. 
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Atom 


WM 








C2(1S) 




1.3336 


3.1031 


0.6502 - 0.0719 


0.4558 - 0.0330 


5.4378 




1.8551 


3.7967 


1.0755 - 0.1525 


0.7492 - 0.0682 


7.2558 




1.9590 


3.9166 


1.1655 - 0.1719 


0.8110 - 0.0765 


7.6038 




2.3937 


4.3693 


1.5575 - 0.2635 


1.0790-0.1150 


9.0209 


ttD 


2.5608 


4.5253 


1.7137- 0.3032 


1.1853 - 0.1315 


9.5504 



To calculate the complete corrections to any quan- 
tity one indeed has to take into account relativistic cor- 
rections and corrections to the operators, which are not 
universal. Let us briefly discuss possible applications. 

Let us consider specifically two systems, namely, the 
muonic hydrogen and pionium. 

(i) For the pionium a part of the corrections is 
known and progress in experiment requires im- 
provement of theory and calculation of remaining 
terms, which are in particular presented in this paper. 

(ii) Another application of interest is due to the hy- 
perfine splitting in muonic hydrogen. In this case the 
complete result includes various nuclear-structure depen- 
dent effects and the corrections are rather irrelevant 
for comparison with ex pec ted experimental data both for 
the Is lll!| and the 2s [l2| hyperfine intervals. However, 
if both experiments will deliver accurate results, one can 
consider a specific difference Q 



AE21 = 8 X Ems{2s) - EmsM 



(5) 



and for this difference the calculation of the terms are 
relevant. The contribution induced by the correction to 
the wave function is only a part of the complete result 
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FIG. 4: The second-order vacuum polarization corrections to 
the hyperfine splitting in muonic hydrogen. The wavy line 
stands for the hyperfine interaction. 

Performing a calculation of third terms depicted in 
Fig. [3Ji and Fig. UJ; as a test, we also calculated a contri- 
bution to the Lamb shift (see Fig. [5]), which was previ- 
ously calculated for muonic hydrogen in (l3 | . 
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FIG. 5: The a m correction to the Lamb shift in muonic 
hydrogen: the only contribution of the third order of non- 
relativistic perturbation theory (see Q). 

The result [3] was for the correction to the Lamb shift 
(i.e. for a splitting of the 2s and 2p states). This value is 
of a particular interest because of the PSI experiment 
p^ . We have also calculated the same quantity and 
found 



AE. 



(Fig.[5l 

2p-2s 



(this work) 



(-7.3861.10-^ + 0.3511.10" 



4 



-(-0.002 5412+0.0013661) 



= 0.0011681^TOr 



(9) 

where the first parenthesis is for the 2p contribution, 
while the second is for the 2s one; each parenthesis con- 
sists of a main term and a subtraction term as introduced 
in ©. 

Our resuh disagrees with the result published in (l3 |. 



the hyperfine structure in muonic hydrogen for the same 
states and corrected the result on the a^m^c^ con- 
tribution to Lamb shift in muonic hydrogen. Details of 
our calculations of the second-order vacuum polarization 
effects for the wave function at origin in various atoms 
and the hyperfine splitting and the Lamb shift in muonic 
hydrogen are in preparation and will be published else- 
where. 



(Fig. m 

2p-2s 



(Ref. [l|) = 0.002535(l)^m^c2 



(10) 



After the work was finished we contacted the authors 
of [l3| . As the results of the communications it has been 
agreed that our result for the main terms for both states 
(2s and 2p) confirm calculations in [l3|, while the sub- 
traction term was missing there. After their result was 
corrected for subtraction, it agrees with ours (see our 
eprint [l^ for detail. 

Concluding, we calculated non-relativistic corrections 
in the relative order to the wave function at origin in 
muonic and exotic atoms for the Is and 2s states, pre- 
sented a result on the non-relativistic correction to 
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